Experimental Measurement of the Persistence Exponent of the Planar Ising Model 
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Using a twisted nematic liquid crystal system exhibiting planar Ising model dynamics, we have 
measured the scaling exponent 6 which characterizes the time evolution, p(t) ~ t~~ e , of the probability 
p(t) that the local order parameter has not switched its state by the time t. For 0.4 seconds to 200 
seconds following the phase quench, the system exhibits scaling behavior and, measured over this 
interval, 6 — 0.19 ± 0.031, in good agreement with theoretical analysis and numerical simulations. 
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There has been a recent surge of interest in determining 
the so-called nontrivial "persistence" exponent 9 jl] Q 
associated with the dynamics of a phase ordering system 
following a quench from the high-temperature phase to 
zero temperature. This exponent describes the asymp- 
totic power-law decay, p(t) ~ t~ e , of the probability that 
the local order parameter (f>(x,t) has not changed sign 
up to time t after the quench. For the Ising model, p(t) 
is simply the fraction of spins that have not nipped up 
to time t, or, equivalently, the probability that no in- 
terface has ever crossed a given spin. While there have 
been considerable theoretical and numerical efforts di- 
rected towards calculating 9, so far there has been no 
experimental measurement of this exponent in spin-like 
systems. A kind of persistence exponent was first intro- 
duced and measured in a breath figure experiment [jjj, 
and recently for soap bubbles |L2| (in this latter example 
9 takes a rather trivial value, which can be inferred on 
simple physical grounds fl^] ). In this Letter, we present 
the first experimental measurement of the nontrivial per- 
sistence exponent for an Ising-like system, that we find 
in good agreement with a recent theoretical calculation. 

This power-law decay of persistence is quite ubiquitous 
and not just restricted to the phase ordering dynamics of 
the Ising or Potts models. For instance, a similar ques- 
tion arises in the study of Gaussian processes. Consid- 
ering such a process m(t), the calculation of the proba- 
bility that this Gaussian walker never crosses the origin 
(or changes sign) is a particularly difficult problem when 
this process is not Markovian [jl0],[ll]]. For example, even 
for the simple diffusion equation, d t <j) — V 2 </>, the proba- 
bility that the Gaussian variable 4>(x, t) does not change 
sign up to time t, decays algebraically with t with a non- 
trivial exponent 9 pi]] . This exponent 9 also appears in 
other contexts such as reaction-diffusion systems 13 15| ] 
and driven diffusive systems |16|. For a quench to the 
critical point of spin systems, the persistence exponent 
9 associated with the total magnetization (another non- 
Markovian Gaussian variable) has recently been argued 
to be a new non-equilibrium critical exponent fli~7| ]. 

For the purpose of the experimental results presented 
below, we will restrict ourselves to the phase ordering 



dynamics of the Ising model at T = 0. The Ising model 
at T = has two ordered states: either all spins are up 
or they are all down. Following a quench from the high- 
temperature initial state where all spins are randomly 
up or down, the system tries to order locally. However, 
since the symmetry between the two ordered phases is 
not broken by the quench process, the two ground states 
compete with each other. As a result, domains of both 
phases form, grow with time, and the system exhibits 
coarsening. At late stages of coarsening, the morphology 
of the growing domains seems to exhibit a self-similar 
structure in time once all length scales are rescaled by 
L(t), where L{t) represents the typical linear dimension 
of a growing domain. The scaling hypothesis also pre- 
dicts that the equal-time two-point correlation function 
has a simple scaling form, (cf>(0,t)(f>(r,t)} ~ g[r/L(t)], 
where the calculation of q{x ) has been the subject of in- 
tense theoretical efforts [|18|. Also, the length scale L(t) 
grows algebraically with time, L(t) ~ t 1 / 2 . This growth 
law can be simply understood by noting that the motion 
of the domain walls between the opposite phases is purely 
curvature driven []l8| . 

The equal-time correlation function, however, only 
gives information about the spatial structure of the sys- 
tem at a given time. To obtain information about the 
temporal evolution, one useful measure is the unequal- 
time autocorrelation function, C(t,t') = (<fi(r, t')0(r, i)). 
For t ^> t', this autocorrelation decays as C(t, t') ~ 
[L{t)/ L(t')]~ x , where A was argued to be a new universal 
nontrivial non-equilibrium exponent Jl9|] . For the d = 2 
Ising model, Fisher and Huse give a heuristic argument 
leading to A = 5/4 |l9| , in good agreement with their nu- 
merical simulations. The experimental measurement of 
A was carried out for the first time by Mason et al. [Q , 
using a twisted nematic liquid crystal sample bounded 
by two glass plates (as will be explained in detail below). 
The experimental value of A in d = 2 was in very good 
agreement with the theoretical prediction |Q . 

The important question, however, is whether this auto- 
correlation is sufficient to characterize the full temporal 
evolution of the system. The answer is presumably no, 
since the stochastic evolution process of a given spin is 
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not Gaussian and not Markovian. What is the minimal 
set of quantities needed to specify the full distribution 
of this complex spatial and temporal structure? It is ex- 
tremely difficult to answer this general question as the 
evolution of the order parameter field in general satisfies 
a nonlinear partial differential equation. In the absence 
of such information, one therefore looks for simple, easily 
measurable quantities that still give important informa- 
tion about the history of the evolution process. One such 
simple and natural quantity is the "persistence", i.e., the 
probability p(t) ~ t~ e ~ [L(t)]~ that a given spin does 
not flip up to time t. 

This fraction of unflipped spins p(t), though relatively 
simple to measure in numerical simulations P,|7|P,|l0| , 
proves extremely hard to compute theoretically. The 
temporal evolution of an individual spin is a non- 
Gaussian, non-Markovian process, and naturally, any his- 
tory dependent nonlocal quantity is very hard to compute 
for such processes. Indeed, this persistence probability is 
not simply a two-point correlation function of the initial 
and final times, but involves all intermediate times as the 
process is, in general, non-Markovian. Nevertheless, for 
the d = 1 Ising model with Glauber dynamics at T = 0, 
this nontrivial exponent 9 was recently computed exactly 
by Derrida et al. who found 9 — 3/8. Unfortunately, 
their method cannot be extended to higher experimen- 
tally relevant dimensions. For the d — 2 spin flip dynam- 
ics of the Ising model, only numerical estimates of the 
exponent 9 w 0.22 |^|,^J|Jl^] have been available to date. 

An approximate analytical method has recently been 
developed |l0| to compute 9 in any dimension for the de- 
terministic nonlinear Landau-Ginzburg equation, which 
describes the T = dynamics of an isotropic spin system 
in the continuum. It involves the following few steps. 
First, one assumes that the late-time dynamics of the 
discrete Ising spins are correctly described by the con- 
tinuum noiseless model-A equation for the nonconserved 
order parameter field </>(r, t). In fact, for the experimental 
system discussed below, this continuum equation of mo- 
tion for the field is certainly more appropriate [ plf . Next, 
one makes a nonlinear transformation (j)(r,t) = cr(m,i), 
where m(r, t) is an auxiliary field locally representing the 
distance to the nearest domain wall |^2|,^3| and a is the 
equilibrium profile of a domain wall, a has a sigmoid 
shape, which, in practice can be treated as the sign func- 
tion. One then assumes that m(r,i) is a true Gaussian 
field (an approximation becoming exact in the large d 
limit |2^,^j|) which allows the calculation of the two- 
point correlators of both m and </> self-consistently. This 
method was originally introduced by Mazenko p2| , p3| |, 
and proves very successful in calculating the exponent A 
or the equal-time correlator g(x). For example, in d = 2, 
this approximation yields A w 1.289 to be compared with 
1.25 ± 0.01 obtained in simulations [Q. Once the above 
approximation has been made, the probability of not flip- 
ping a spin is the same as the probability that the Gaus- 



sian process m, at a given point in space, does not change 
sign up to time t, as the equilibrium interface profile is 
an odd function (0 = er(m) ~ sign(m)). Still, the com- 
putation of the probability of no zero crossing for the 
Gaussian process remains extremely hard due to its non- 
Markovian nature. 

However, the problem of computing 9 can be mapped 
onto a quantum mechanics problem [ p)[ of sorts. 9 is then 
found to be the difference between two ground-state ener- 
gies: that of a quantum-like problem with a hard wall at 
the origin (this constraint results from the condition that 
the Gaussian walker has to remain on, say, the positive 
axis), and the other for the same system but without the 
wall. In these terms, the quantum problem associated 
with a Markovian process is simply an harmonic oscil- 
lator with frequency A/2, where A is again the autocor- 
relation exponent. One can then develop standard vari- 
ational and perturbative methods to compute the expo- 
nent 9 for a process close enough to a Markovian process. 
In d = 1, our approximate theory yields 9 « 0.35 Jn| as 
compared with the exact value 6> = 3/8 Q. For d = 2, 
we found 9 « 0.19 G3|. This last result is also consistent 
with the estimate 9 « 0.186 obtained for d — 2 in the 
asymptotically exact d — ► oo limit of Mazenko theory us- 
ing the independent intervals approximation for the zeros 
of the process m(t) pl] |, and in recent direct simulations 
of the time-dependent Ginzburg-Landau (TDGL) equa- 
tion (model-A) |2lJ]. Note that if the process m(t) were 
Markovian, one could show [|l0| that the simple equality 
9 = A/2 w 0.63 holds. The discrepancy with the actual 
value is an indication of strong memory effects in this 
system. 

For comparison with theory, liquid crystals afford the 
possibility of designing experiments for which 9 can be 
easily measured. Indeed, a number of liquid crystal 
systems, possessing a variety of dimensionalities and 
types of order parameters, have been used to mea- 
sure various scaling exponents characterizing the coars- 
ening process. Such measurements have included the 
increase in the characteristic scaling length, L(t), in 
three-dimensional J|4|-^7j and two-dimensional |28|-|3"2j| 
systems. The twisted nematic system, exhibiting Ising- 
like behavior and employed here for the measurement of 
9, has been previously used by Orihara and co-workers 
p8| , p9| to study the behavior of the scaling length L{t) 
and, in addition, by Mason et al. [^(J to measure A. These 
studies used a liquid crystal sample placed between two 
glass plates whose surfaces had been treated so that, after 
a thermal quench from the isotropic phase to the nematic 
phase, the liquid crystal organized itself into domains in 
which the director was forced to twist clockwise or coun- 
terclockwise by it/ 2 in going from the surface of one plate 
to the other. The boundary between two domains of op- 
posite twist consists of a twist disclination line. Regions 
of opposite twist correspond to Ising model domains in 
which the spins all point up or all point down. The sys- 
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tern relaxes viscously, driven by the line tension of the 
disclination lines. 

The experimental apparatus used for the A mea- 
surements has been described in detail by us previ- 
ously and the data collected for this report was 
extracted from the video tapes obtained in these ear- 
lier Ising system studies. Briefly, the sample cell con- 
sists of a layer of liquid crystal, 20^m thick, sand- 
wiched between two parallel glass microscope slides. 
The liquid crystal used was trans-(trans)-4-methoxy-4'- 
n-pentyl-l,r-bicyclohexyl (Merck, CCH-501, or equiv- 
alently ZLI-3005), with an isotropic-to-nematic phase 
transition of 37 °C, The molecular orientation at each 
slide's surface was forced to be in the plane of the surface 
and oriented unidirectionally. The surfaces of the glass 
plates were prepared by first dipping them into a 0.1% by 
weight solution of polyvinyl alcohol and then buffing the 
surface in one direction with a soft cloth. The two slides 
were mounted orthogonally to each other, thus defining 
a square region of liquid crystal, about 2.5 cm on a side. 
The orthogonal orientation of the plates with respect to 
one another forces the nematic liquid crystal to phase 
separate into two domains, as the molecular orientation 
undergoes either a clockwise or counterclockwise twist of 
7r/2 radians as one goes from one plate to the other. The 
circular dichroic effect was used to make the domains 
of opposite twist visible. A square region, 1.5 mm on 
a side, was observed using a Nikon E Plan x4 objective 
and recording the images with a high-speed color record- 
ing system (NAC, HSV-400). Images were recorded every 
5 ms, each frame labeled at the top with the run number 
and time since the recording began. 

For a well percolated system undergoing self-similar 
coarsening, the ratio of the area occupied by domains of 
one particular twist to the total area as a function of time 
should remain constant and have the value 0.5. This ra- 
tio, averaged over the 15 runs discussed here, was 0.375 
at the time of the quench, reaching within only 0.2 sec 
the value of 0.576 ± 0.016, at which it remained for the 
200 second time interval over which data was recorded. 
There are several possible reasons for the deviation of 
this ratio from the value of 0.5. The region between do- 
mains over which the director field is distorted from the 
pure clockwise or counterclockwise twist is comparable 
to the thickness of the cell, 20 microns. This makes the 
optical image of the boundary between domains some- 
what ill-defined. Hence, optical effects could cause us 
to preferentially choose one domain over the other when 
deciding whether a given point belongs to a region of 
clockwise or counterclockwise twist. This would be par- 
ticularly serious at early times and probably accounts for 
the low value (0.375) of the ratio measured at early times. 
The ratio at late times is 15% higher than the expected 
value of 0.5 and may be due to a small bias in the do- 
main nucleation, domain growth favoring domains of one 
particular twist, and, at very late times, small sample 



statistics. For example, if alignment of the upper and 
lower plate deviate from perfect orthogonality one twist 
direction will be favored over the other |2i| . 
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FIG. 1. Probability p(t) versus time that an individual 
point in the system has not switched its state at least once 
by time t. The solid line is a least-squares fit to the data and 
has a slope of 0.19. 

We measured, for each of 40 points (in an ordered, 
rectangular array) in each of 15 runs |33|, the time t after 
the quench when the molecular orientation first switched 
from one twist direction to the other. The contrast was 
sufficient to enable us to distinguish the orientation at a 
particular point immediately after the phase transition. 
Fig. 1 is a log-log plot showing the scaling of the prob- 
ability pit) an individual point has not switched its ori- 
entation by time t. A least-squares fit to the data in the 
range 0.4 to 200 seconds gives a slope, 8 = 0.190 ±0.031, 
shown by the solid line. The initial approach to scaling 
is most likely due to the initial nonzero size of the scale 
length L at the time of the quench and the identification 
bias resulting from the diffuse nature of the optical image 
of the domain boundaries. This initial scale length is on 
the order of the thickness of the cell (20/mi), resulting in 
a minimum domain size at the time of the quench. Also 
noticeable at late times is a "tail" , of reduced slope, in 
p(t). The tail at late times could be due to several differ- 
ent effects: pinning of domain walls in the sample imme- 
diately outside of the region observed or nonorthogonal 
orientation of the glass plates (giving rise to a preferred 
molecular twist). The wiggles in the data are finite sam- 
ple size effects reflecting the history of how macroscopic 
regions changed their orientation in particular runs. 

In conclusion, our measured value of the persistence 
exponent 8 = 0.190 ± 0.031 is in good agreement with 
the theoretical estimate 8 w 0.19 for an isotropic bidi- 
mensional spin system. 
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